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Abstract. - We have computed the time-dependent susceptibility for the finite-size mean-field 
Random Orthogonal model (ROM). We find that for temperatures above the mode-coupling 
temperature the imaginary part of the susceptibility x' ' W) obeys the scaling forms proposed 
for glass-forming liquids. Furthermore, as the temperature is lowered the peak frequency of 
x" decreases following a Vogel-Fulcher law with a critical temperature remarkably close to the 
known critical temperature T c where the configurational entropy vanishes. 



A model Hamiltonian or an effective Lagrangian capable of describing relaxation processes 
in supercooled liquids and structural glasses is difficult to obtain. However, starting with 
the work of Kirkpatrick, Thirumalai and Wolynes in the late 80's, it is now clear that 
there is a close analogy between some mean- field spin-glass models and structural glasses [3] . 
The former thus provide a set of microscopical models where glassy dynamics can be studied 
analytically. The basic simplification occurring in mean-field models is that after averaging 
over the disorder and making the number of spin very large (N — > oo) one is left with a closed 
set of equations for the two-time correlation and response functions which, above a critical 
temperature, are basically equivalent to the schematic mode coupling equations introduced 
by Leutheusser, Gotze and others |3] as a model for the ideal glass transition. 

In mean-field models the barrier separating different ergodic components diverges in the 
mean-field limit, hence at the critical temperature To a real ergodic to non-ergodic transition 
takes place with diverging relaxation times since barriers cannot be overcome, in agreement 
with mode-coupling theories (MCT) where To coincides with the MCT critical temperature 
Tmct- In what follows we will denote To by Tmct- In real systems Tmct still denotes the 
separation temperature for activated hopping dominated dynamics [5], however, barriers are 
of finite height and the glass transition appears at T g < Tmct where the typical activation 
time over barriers is of the same order of the observation time. To go beyond mean-field it is 
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necessary to include activated processes, a very difficult task since it implies the knowledge of 
excitations involved in the dynamics. 

Recently it has been shown how activated processes in mean-field models could be included 
just keeping N finite [$]], and hence giving support to the scenario of the fragile glass transition 
developed from spin-glass models. This is not a trivial assumption since it is not a priori clear 
why excitations in mean-field spin glass models should have similar properties to those of 
supercooled liquids. This, for example, seems to be the case for the Random Orthogonal 
model (ROM) [7J, but not for the mean-field Potts-Glass model as indicated by recent studies 
of its finite-size version [5] . 

In a series of works the spectral properties of the primary or a relaxation in supercooled 
liquids has been largely studied by means of dielectric spectroscopy |^l ll(JlliTlll2H13j . The 
main conclusion is a there is a three-parameter scaling function that allows to collapse the 
imaginary part e"(v) of the dielectric susceptibility e(v) at different temperatures and for 
several glass-forming liquids onto a master curve. The master plot is able to reproduce the 
e"(i/) data around the relaxation peak v p and also at higher frequencies. The bad collapse in 
the low frequency part has been object of some debates |1()II12| . however there is no dispute 
above v p . The frequency v v has a very strong temperature dependence, commonly fitted by 
a Vogel-Fulcher form log 10 (^ p ) = log 10 (^ ) — A/(T — T ), where T is close to the Kauzmann 
temperature |14| where the configurational entropy vanishes |15j . 

In this contribution we compare the frequency-domain analysis of the finite-size ROM 
with the above scenario. The main motivation for this study was to make a stringent test 
on the ROM as a possible model for the fragile-glass scenario. The range of temperatures 
we explore in this paper are all above the mode-coupling transition Tmct, therefore we do 
not expect to find diverging timescales in the large N limit. In the light of the experimental 
results 0E3^nElEj, that also explore the region above Tmct, we have decided to focus 
our attention in this regime in the framework of finite-size mean-field models. Here, we do 
not address the behavior below Tmct where a N dependence of the frequency response is 
expected. 

The ROM [7J is defined by the Hamiltonian 



where cr, = ±1 are N Ising spin variables, and is a N x N random symmetric orthogonal 
matrix with J a — 0. This model was proposed in [7] as a random extension of some non- 
disordered models describing low autocorrelation binary sequences. Numerical simulations 
are performed using the Monte Carlo method with the Glauber algorithm for temperatures 
in the range 0.6 up to 2.0. For N — ► oo and h = this model has a dynamical transition at 
Tmct — 0.536, and a static transition at T c = 0.256... [7J. To study the frequency response 
we considered a time-dependent field of the form h(t) — ho cos(2 7r i/t), where the time is 
measured in MC steps and ho = 0.2 small enough to be within the linear response regime. In 
our simulations the typical range of v was 10 -6 - 10 _1 . For each frequency we have measured 
the complex susceptibility x{ v ) — x'{ v ) + as 
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The number of MC steps M after equilibration was 100 for the largest v and up to 10 7 
for the shortest v. We show here results for N — 300, the finite-size dependence will be 
discussed below. The choice of N = 300 is a good compromise between small sample-to- 
sample fluctuations and small barriers height. 

I 1 I 1 I 1 I 1 I I 




Fig. 1 - The real and imaginary part of the complex susceptibility \' an d x" as function of v for the 
ROM with N — 300 at different temperatures. Upper panel top to bottom, lower panel left to right, 
T = 0.7, 0.9, 1.1 and 1.3. 

Figure ^ shows the real and imaginary parts of the susceptibility over the available range 
of frequency. Not all temperatures are reported for a better drawing. The relaxation peak in 
the imaginary part can be fitted with a log-normal form |16) : 

X» = 4| exp [-(log 10 i/ - log 10 i/ p ) 2 /E 2 ] (4) 

where v v is the frequency of the peak, £ the width, and Ax — Xo ^ x'ooi where x'o an d x'<x> 
are, respectively, the low and high frequency limit of x'(f). 

As the temperature is lowered the peak frequency v v decreases, and the width £ broadens. 
The behavior of v p is consistent with the Vogel-Fulcher law exp[— A/ (T — To)] The fit of 
the frequency peak v p for the ROM with the Vogel-Fulcher formula is rather good, see Fig. 
and gives A = 0.89 ± 0.06 and In v = 0.64 ± 0.02 T = 0.28 ± 0.02, a value in agreement with 
the critical value T c = 0.256.... We note, however, that data can also be fitted using different 
expressions such as exp[— A/T 2 } or the Adam-Gibbs formula exp[— A/ (TS C (T))] where S c is 
the configurational entropy. In particular we note that the data can be also fitted with the 
formula (CR) v v = v exp[— A(3 c g (T)/T] where (3 c s(T) = 9S , c (ei s )/<9ei s | eia=ei3 (T) derived from 
a cooperative scenario of relaxation |17j . see Fig. El The discrepancy at high temperature is 
probably due to a poor numerical estimation of the configurational entropy, indeed a similar 
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deviation is found using the Adam-Gibbs formula (not reported). We note that the CR 
formula predicts a crossover form fragile to strong behaviors as the temperature is lowered. 
However differences among all the above expressions can be appreciated only for very low 
values of v v which are out of our measurements range. 




Fig. 2 - log 10 ^ p as function of T for the ROM with N = 300. The full line is the Vogel-Fulcher law 
u p = v exp[-A/(T - T )] while the dashed line is the formula v p = u exp[-A(3 cB (T)/T] of Ref. [T7] 
with /3 e ff(T) = 9S'c(ei s )/9eis|e is =e is (T) evaluated using the results of Refs. [Sj. 

The analysis of the response for glass-former liquids reveals three power laws for \" .13 : 
t^" ~ v m for v < Up, x" ~ v ~P just above v v and x" ~ 1, f° r larger values v ^> v p . The use 
of Monte Carlo dynamics prevents us from resolving the last behavior due to the discreteness 
of the time step, nevertheless the first two regimes are clearly seen, as indicated in Figure [3] 
At higher temperatures m = (3 = 1 and the relaxation is Debye-like. As the temperature is 
lowered the value of j3 decreases below 1. It is known that for glass- former liquids j3 and a are 
related by (er + + 1) = 7, where 7 is a constant Furthermore a varies linearly with 
temperature: a = B(T - T a ) with T CT ~ T This implies that = B'(T - T ) + (1 - 7 )/ 7 . 
Inserting into this formula the values of (3 obtained for the ROM at various temperatures 
and the value of To computed from v p we find 7 = 0.72 ± 0.02 the same value found for real 
liquids EMEU- 

Analyticity of x( l/ ) an d linearity of absorption at asymptotically low frequencies implies 
that x"{ v ) v f° r v ^ v v CHI- For the ROM with N = 300 we find m ~ 1 for temperatures 
down to about 0.8. Below T = 0.8 significant deviations with m < 1 are observed. The 
temperature where m — 1 breaks down decreases as N increases. Similar deviations have 
been observed in data from glass-forming liquids and generated some controversy [91 1101112] 
on the reliability of the scaling form proposed by Dixon et al. [H]- Many liquids also posses 
secondary relaxations which overlap the primary response broadening the peak and leading to 
deviation from linearity ^2]. In the case of ROM these secondary relaxations are related to 
the fact that the barriers separating the low states sampled as the temperature is decreased 
toward Tmct are not well separated for not too large N. Indeed studies of mean-field spin- 
glass models for the structural glass transition shows that in the thermodynamic limit there 
is no gap between saddles separating local minima with energy above the threshold energy 
associated with the dynamical transition JSj- This is a situation more reminiscent of spin- 
glasses rather than glasses. Indeed both experimental an d numerical simulations |19U2U| 
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Fig. 3 - x' i u ) f° r the ROM with TV — 300 for temperatures T = 0.7 and 1.1. The exponents are 
m = 1, P = 1. for T = 1.1 and m = .96, j3 = .71 for T = 0.7. Inset: *"(» at temperature T = 0.6 for 
TV = 64 (filled circles) and 300 (empty circles). The lines have slope m = 0.3 for TV = 64 and m = 0.5 
for TV = 300. The increase of m toward 1 as TV grows is clearly seen. 



show a broader shape of x"{ v ) near the peak. 

This scenario is supported by a finite-size scaling analysis of the ROM. Indeed we find 
that for a fixed temperature while (3 is independent on TV, the value of m, when less than 
1, increases toward 1 as TV is increased. In the inset of Figure |3| we show x"{ v ) at T = 0.6 
for systems of size TV = 64 and TV = 300. The increase of m going from 64 to 300 spins is 
well evident. We now finally address the goodness of the Nagel scaling. Dixon et al 9 have 
shown that all data for the dielectric susceptibility e{v) of different glass-forming liquids and 
temperatures can be collapsed onto a single master curve by plotting (l/io)log 10 (e"(z/)v p /VAe) 
versus (l/w){X/w+X)\og w {y/v p ), where Ae = e'(0)— e'(oo) is the relaxation strength measured 
by the step height of the real part of e(y), v v is the peak frequency, and w is the half- 
maximum width of the e"{v) peak normalized with the corresponding width of the Debye 
peak: Wd — 1-14 decades. 

FigureElshows the Nagel plot for the ROM with TV = 300. The data are from temperatures 
ranging from 0.6 up to 2.0. For each curve the parameters Ax, w = 2v / 2£/u; ( i and v v have 
been obtained from the log-normal fit of x" using J3J- We see that while the collapse for 
v > v p is good for all temperatures, for v < v v only data with m — 1 do collapse. As noted 
in Refs. the optimization of the three parameters through the fitting is essential to have 
a good collapse of data. For comparison in the inset we report the Nagel plot for the same 
data done using the parameters read directly from the plot of Figure ^ without any fitting. 
The improvement through the fitting is evident. 

In conclusion, we have shown that the primary relaxation in the finite size mean-field 
ROM obeys the scaling form typical of glass-forming liquids. Furthermore, the frequency 
peak of the imaginary part of the complex susceptibility follows the Vogel-Fulcher law with 
critical temperature T = 0.28 ± 0.02 very close to the critical temperature T c = 0.256..., 
the Kauzmann temperature of the model. All system sizes studied (up to TV = 300) lead 
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w"'(l+w"') log (v/v 



Fig. 4 - The Nagel's scaling for the ROM with TV = 300 and temperatures T = 2.0, 1.3, 1.2, 1.1, 1.0, 
0.9, 0.8, 0.7 and 0.6. See text for details. 



to this value for Tq. Because we used Monte Carlo dynamics there is a maximum value 
for the frequency ~ N determined by the discreteness of the elementary time step. As a 
consequence we are not able to resolve the second scaling behaviors of the imaginary part of 
the susceptibility x"{ v ) ~ V<T f° r v ^ v p . Nevertheless by assuming that the exponent a 
vanishes linearly at To, we obtained for the constant 7 relating the exponents /3 and a the 
value 7 = 0.72 ± 0.02, the same found for real glass-forming liquids |11II13| . Again this value 
does not depend on the system size for the sizes we studied. Changing N only modifies the 
temperature below which secondary relaxations related to barrier crossing breaks the Nagel 
scaling. Overall, the present results show that finite-size mean-field spin glasses capture the 
cooperative effects responsible for the relaxational processes observed in glass forming liquids 
when approaching the mode-coupling temperature. Moreover, to our knowledge, this is the 
first time that finite-size mean-field spin glasses give a quantitative prediction in agreement 
with experimental data. The extension of this analysis to the region below Tmct where strong 
finite N effects arc to be observed remains an interesting open problem. 
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